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Abstract 

We describe the quantum dynamics of the Hubbard model at semi-classical 
level, by implementing the Time-Dependent Variational Principle (TDVP) 
procedure on appropriate macroscopic wavefunctions constructed in terms of 
su(2)-coherent states. Within the TDVP procedure, such states turn out to 
include a time-dependent quantum phase, part of which can be recognized 
as Berry's phase. We derive two new semi-classical model Hamiltonians for 
describing the dynamics in the paramagnetic, superconducting, antiferromag- 
netic and charge density wave phases and solve the corresponding canonical 
equations of motion in various cases. Noticeably, a vortex-like ground state 
phase dynamics is found to take place for U > away from half filling. More- 
over, it appears that an oscillatory-like ground state dynamics survives at the 
Fermi surface at half-filling for any U. The low-energy dynamics is also exactly 
solved by separating fast and slow variables. The role of the time-dependent 

phase is shown to be particularly interesting in the ordered phases. 
71.10.Fd,03.65.s,03.65.B,47.32.Cc 
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I. INTRODUCTION 



Interest in strongly-correlated itinerant electron systems has been constantly growing 
in the last three decades. Especially since the discovery (almost ten years ago) of high-T c 
superconductors, an enormous amount of work has been devoted to such systems, aimed both 
to investigate their macroscopic thermodynamical properties via experimental mesurments, 
and to disclose -by employing the standard methods of statistical mechanics- what type of 
macroscopic collective order is responsable for the frictionless regime. 

Nevertheless, due to the high number of variables naturally involved in the models pro- 
posed for investigating these many-electron systems, and probably to the background of the 
community of physicists who first considered these models, to the best of our knowledge 
very little effort has been made in order to investigate their dynamical behavior. On the 
other hand, this type of analysis is known to lead to interesting properties of superfluidity 
when applied for instance to the BCS HamiltonianS. 

Two circumstances, at least, prompt to attempt the dynamical approach and to carefully 
consider its possible developments. First of all, standard theoretical techniques such as the 
Time-Dependent Variational Principle (TDVP) procedure and its path-integral version, the 
Stationary Phase Approximation method, have been remarkably developed in the recent 
years, by exploiting the notion of Generalized Coherent State (GCS)i and the spectrum 
generating (or dynamical) algebra method. At the formal level, such group-theoretical tools 
have greatly simplified and provided of a systematic character the TDVP procedure, which 
essentially consists in reducing the system quantum dynamics to a semiclassical hamiltonian 
form. The procedure, formulated, for example, in the form of refs.i and!, was first introduced 
for studying the low-lying collective states in nuclei, but it is easily extended also to any 
systems endowed with a large number of degrees of freedom. 

Furthermore the special role assigned to the quantum phase of the macroscopic trial 
wavefunction involved within the TDVP framework makes the procedure even more attrac- 
tive. Such quantum phase, in fact, is structured so as to have a memory of the entire 
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dynamical behavior. It is thus natural expecting some kind of correlation between the type 
of microscopic order which possibly characterizes the medium and the phase time-behavior. 
Such kind of effects have been for instance investigated inl where the study of quantum 
dynamics of solitons in ferromagnets was shown to have remarkable consequences on Berry 
phase behavior. 

This aspect, in turn, directly leads to the second circumstance which motivates our inter- 
est for the dynamical viewpoint. At low temperature, one can reasonably think of performing 
current measurements or superconductive quantum interference measurments based on an 
experimental devices similar to the ones employed to investigate the Josephson effect^. It 
is essential recalling that, in that case, the time-dependence of the order paramenter (the 
macroscopic wavefunction of the system) has a crucial role in allowing for the detection of 
microscopic phenomena occuring in the medium. Although the TDVP approach is able of 
taking into account large number of dynamical degrees of freedom, and is thus fitting to 
describe a strongly-correlated electron gas, one can expect that low excited states actually 
involve a restricted number of dynamical variables. Under such conditions the quantum 
phase could retain a nonrandom character which makes it suitable for experimental measur- 
ments. 

A further quality of the TDVP approach must be still pointed out. Such scheme, in 
fact, involves the construction of new semiclassical Hamiltonians which are obtained in a 
completely independent way with respect to the model Hamiltonians derived by standard 
mean-field techniques of statistical mechanics (for instance Hartree-Fock). The ground state 
of the semiclassical Hamiltonian reproduces, as we shall see, the same results of the Hartree- 
Fock approach from the set of dynamical fixed points. Moreover, as opposite to the mean- 
field cases, here also the excited states at low energy are expected to be a realistic description 
of those of the original model Hamiltonian. This character is related with the fact that TDVP 
Hamiltonians, even though affected by the approximations imposed by the method, generally 
preserve a structure rather faithful to the second quantized Hamiltonian. Comparison with 
statistical mechanics approximate models is thus interesting in any case. 
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In this paper we propose the implementation of an appropriate generalization of TDVP to 
itinerant interacting electron systems. This amounts to applying the TDVP to a trial wave- 
function representing a semiclassical macroscopic state constructed by generalized coherent 
states of the dynamical algebra of the model Hamiltonian&i. From the semiclassical picture 
of the system obtained in this way one can derive canonical equations of motion, and a 
classical description of system's dynamics. A key role within this approach is played by 
a time-dependent phase, which has to be fixed so as to satisfy -at least in average- the 
Schrodinger equation. In@ it was shown that, under appropriate assumptions, the latter is 
nothing but the dynamical plus geometric phases! beyond the adiabatic approximation. 

Here we apply the above method to the Hubbard mo dei, described by the Hamiltonian 

H H ub = ~T Yl c t° c ^ + U Y1 ni >T n U - ^ $^Kt + n u) 

<i,j> cr i i 

= X^ Ck ~ V) 71 ** + U Yl a k,T a ma a U a k+m-l,T , (1) 
k.cr k,l,m 

where, on the first line, c| CT ,Cj i(7 are fermionic creation and annihilation operators 

({cj j(T , Ci >(7 '} = , {cj j(J , ci^,} = 5ij S ffj(T rI, n- hC r = <\ a c- h(T ) on a d- dimensional lattice A (i, j G A, 

o G {it, Jj}) with N sites, and < i, j > stands for nearest neighbors (n.n.) in A. In the second 

line the same Hamiltonian is rewritten in the reciprocal space A, with = e CT -*' k Cj iCr , 

j 

d 

Ck = — 2T^^cosfc r . In (|I]) the first term represents the tight-binding band energy of the 

r=l 

electrons (T being the hopping amplitude), the U term describes their on-site Coulomb in- 
teraction, and fi is the chemical potential, which will allow us to fix the conserved quantity 
M e = ^2i( n i,i + n u)> i- e - the total electron number operator on the lattice. 

Being the GCS of the dynamical algebra of Hamiltonian quite complex to deal with, 
we choose as trial GCS's for constructing the semiclassical macroscopic state the su(2) 
coherent states which are exact for the corresponding Hartree-Fock Hamiltonian. This is 
done for two different cases, namely that describing superconducting (SC) and paramagnetic 
phases, and, at half-filling, that describing antiferromagnetic (AF) and charge-density-wave 
(CDW) phases. The approximate equations of motion we obtain for the full Hamiltonian 
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in the two cases are then solved in some integrable cases, and by approximate methods in 
other interesting limits. This gives rise to a variety of different dynamical behaviors, from 
vortex-like dynamics in the ground state for the U > 0, n / 1 regime, to oscillations of 
the number of electrons around the Fermi surface and possible laser effect at low energy, 
and to single-mode collective frequency dynamics, which should reflect the occurence of 
macroscopic order in the medium. In particular, the time-dipendent phase, which - due to 
its macroscopic nature - can be considered as an observable quantity, exhibit a behavior 
which is shown to be related to the non-vanishing of order parameters, and is evaluated 
exactly in several situations. 

The paper is organized as follows. In section II we review the generalized TDVP approach 
and its connections with the quantum geometric phase. In section III we treat explicitly 
the SC-paramagnetic case, by constructing first the macroscopic trial wave-function and 
the corresponding semiclassical Hamiltonian, and deriving then the canonical equations of 
motion with the time-dependent phase factor. Section IV is devoted to study the fixed 
points of these equations and in particular the ground state solutions and metastable states; 
obtaining for the ground state the Hartree-Fock results as well as a nontrivial vortex-like 
dynamics in the so called paramagnetic phase, and topological excitations for the metastable 
states. In section V we analyze an integrable case, which exhibits collective order and non- 
zero pairing induced by the /c-modes interactions. In section VI we investigate the global 
dynamics by comparing slow with fast degrees of freedom. We show how slow variables 
tend to constitute an autonomous subsystem which drives fast variables dynamics on large 
time-scales. An integrable case where the slow subsystem is reduced to a two-level system 
representing A;-modes close to the Fermi level, is explicitly solved in section VII. In section 
VIII we repeat some of the above analysis for the AF phase at half-filling, finding in particular 
an oscillating behavior at the Fermi surface. The final section is devoted to some conclusions. 
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II. GENERALIZED TDVP METHOD 



Knowledge of the dynamical algebra Q of a given (time-independent) Hamiltonian H 
allows the construction of an over-complete set of states known as generalized coherent 
states, 



|*q >= exp \Yj[e a E a - e* a E- a ] \ |0 > 



where {E a ,E^ a } are the raising and lowering operators in the Cartan representation of Q, 
and |0 > is the highest weight of the representation, defined by E_ a \0 >= for all positive 
a's. The state >, which is the obvious time-dependent generalization of |$ >> 

|$(t) >= exp (y2[6 a {t)E a - e a (t)*E-a] \ |0 > , 

is related to the time evolution of |$o >, described by the state \^{t) >= e — ^ | <^> >, 
through 

\V(t)>=e i£ % l \$(t)> , (2) 

where 

p^-nt + ih dr < $(r)| — |$(r) > , (3) 
Jo v T 

and H =< $(*)|#|$(t) >. 

The time-dependence of the parameters #a,(i)'s is determined by imposing that \^(t) >, 
as given by (§), satisfies the Schrodinger equation. It turns out that this amounts to requiring 
the a (t)'s obey the canonical equation of motion (see below). An alternative parametriza- 
tion for the state |$(t) > can be used, namely 



|$(t) >= AMexp [ z a(t)E a J |0 > 

\a>0 / 



(4) 



M being the normalization factor, such that < $(£)|$(£) >= 1. In this case the canonical 
equations of motion readB 







dU 



&H 



(5) 



where g a $ 



d 2 log AT 

dz*dz(3 



is the metric of the phase space spanned by {z a , z*}, and determines 



its symplectic structure. The metric g a> p indeed determines the explicit form of Poisson's 
brackets, 



where g~ x represents the inverse matrix of g. 

If the algebra Q is the full dynamical algebra of H (i.e. if H 6 Q for any choice of 
the physical parameters) the above procedure is exact. In particular, it gives the exact 
quantum ground state of H as fixed point of the equations (||). Nevertheless, in a many- 
body problem like the one described by the Hamiltonian ([!]) the dimension of the dynamical 
algebra is exponentially growing with N, and infinite in the thermodynamical limit. Even 
though one may still work out the canonical equations of motionll!, their explicit solution 
becomes then quite hard to handle. It is therefore reasonable to inquire to which extent the 
above scheme can be used in the case where the GCS are built in a sub-algebra A C Q for 
which the equations of motion become tractable. In this case of course \^{t) >, as given 
by (0), differs from e~% Ht \$o >> an d the Schrodinger equation is in general not satisfied. 
Nevertheless, the answer given by the generalized TDVP approach!! is that in fact the above 
scheme still holds also in this case, if one just requires that at least the inner product of the 
Schrodinger equation for \^f(t) > with < ^(t)\ vanishes, i.e. 



Notice that now \^(t) > is to be built only with the raising operators E a G A. Hence 
the analogy with the exact solution (where E a e Q) is complete at a formal level, 

but approximate in the results, H being evaluated on a sub-space of the whole dynamical 
algebra Q. Nevertheless, as for the reliability of the method, one should recall that in the 






(7) 
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limit h — > the results obtained within TDVP become exact (to second order in h) , meaning 
that one obtains the classical description of the system dynamics. 

Both in the exact case, and in the TDVP approximation, the role of phase <p(t) given by 
(||) is particularly simple at the fixed points of ([|). In fact in this case \&(t) > is independent 
of t, and the second term at the r.h.s. of (the so called kinetic term) is vanishing. This 
implies that \^f(t) >= exp (— — Ht)\$o > is uniquely determined by the energy pertaining 
to the initial state |$o >• Such behavior is very reminiscent of what is called the dynamical 
phase for a time-dependent Hamiltonian in the adiabatic approximation. More generally, 
by inserting (|3]) into the expression (|2|) for the state \^f(t) >, the latter can be written as 

*(*)> , (8) 

and we recognize a formal analogy between the phase ip(t) and the dynamical plus geometric 
phase in the adiabatic approximation (for a derivation of the quantum phase beyond the 
adiabatic approximation, see for instance ref. [10]). More precisely, in (|8]) we can identify the 
first term in the exponential with the dynamical phase, and the other (kinetic) term as the 
geometric phase obtained by relaxing some of the hypotesys of the adiabatic approximation. 
We recall that the latter is nothing but the so called Berry0 phase. In fact in some simple 
exactly solvable caseB it was shown that -by imposing appropriate quantization condition- 
the phase (§) does coincide with the geometric plus dynamical phases even if evaluated within 
the generalized TDVP approximation scheme. This leads us to expect that the present 
approach, apart from leading to a simplified if approximate description of the dynamics of 
the Hubbard model, may give a precise physical information, i.e. which is the Berry phase 
of the states we are studying. If this is the case, we expect that whenever we shall obtain 
states with the same energy 7i but different phases <f(t), an appropriate physical device 
should be able to observe their interference. 

A further investigation of this relationship for the model discussed here is beyond the 
purpose of the present paper. Here we want just to emphasize that the state \^f(t) >, 
thanks to its phase (f(t), is able in principle to approximate the wave-function of the full 
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\^(t) >= exp 



- / *(- 



-H + ih<$(t)\%-\$(t) >) 
at 



model Hamiltonian, no matter how small is the subalgebra A of Q. However, we expect 
in general that the results will be the more reliable the more A is a reasonable description 
of Q, i.e. the more the Hamiltonian is in a (thermodynamical) phase in which the relevant 
operators are contained in A. It will turn out that precisely (p(t) will measure how far 
the system is from the states generated by A. In particular, we expect that the system is 
correctly described by A whenever tp(t) happens to be linearly increasing with time, in that 
the wave-function \^(t) > which satisfies (^) results to differ by the one we constructed in 
A (|$(t) >) just for an oscillating phase factor. 

III. SEMICLASSICAL EQUATIONS OF MOTION IN SC PHASE 

The model described by (0), has been intensively studied in the literatureEl. However, as 
for most many-electron problems, it is quite difficult to obtain rigorous results (for a recent 
review, seeEl). In particular, only the one-dimensional zero-temperature energy is known 
exactlyS Therefore, different approximation schemes have been used in order to deal with 
(|I|). Among them, a standard one is the Hartree-Fock decoupling procedure, which amounts 
in approximating the interaction term by reducing it to a sum of its bilinear parts, weighted 
by coefficients which have to be fixed self-consistently. There are of course different ways of 
decoupling the interaction term, depending on the phase which has to be investigated^!-! 
The dimension of the dynamical algebra of the resulting "decoupled" Hamiltonian turns 
out to be greatly reduced, and to all effects one is lead to deal with a subalgebra of the 
spectrum generating algebra of the original Hamiltonian. It is our purpose to construct the 
GCS involved by TDVP scheme in these subalgebras. 

In order to identify the subalgebra, we explicitly need the reduced Hamiltonian. Neglect- 
ing for the moment the possibility of an AF or a CDW phases (which will be investigated 
in a later section) and of a ferromagnetic phase (which is not to be expected at low U), it 
turns out that such Hamiltonian coincides with the linearized Hamiltonian H\ = H^ c \ 
where 
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H ( k sc) = (e k - ix)n k + U [-n k + (Aa[a T _ fc + h.c)\ - U(—N + |A| 2 )iV . 

Here, as customary, k stays for the multi-index (k, a) (k = (k, f) , —A; = (— k, |)), and 
n k = a\a k . Moreover n is the average electron number per site, n = — < Af e >, and 
A = — < a~k a k > is the average pairing per site, where < • > denotes the expectation 



t„t 



2 



k 

value of operator • over appropriate states. 

It is important noticing that h[ sc \ contrary to H Hub , for any A ^ does not commute 
with the electron-number operator per-site J\f e . This is justified by observing that is a 
faithful approximation of Hu u b in an ordered phase which does not conserve such quantity 
(namely, the superconducting phase). On the other hand for A = Af e is still conserved, 
and H[ SC ^ describes in that case the system in a paramagnetic phase, which is known to be 
the case for the ground state, at least for low positive U, away from half filling (i.e. 1). 

The Hamiltonians H^'s have the property that [H^, H^} = 0, and hence can be di- 
agonalized simultaneously. More precisely, h[ sc ^ turns out to be an element of the dynamical 
algebra A sc = (B k A k c \ where is the local su(2) generated by 



4 SC) = {4 k) = \{n k + n- k - 1), Jf = a[al k , J W = a_ fe a fc } 



(9) 



Any eigenstate of the Hamiltonian can then be expressed as superposition of the GCSs 
of A sc , namely 



\r) >= II fc (l + r) k r) k ) 2 exp (r] k j| fe) j |0 



>, 



where the coefficients i] k e C parametrize the (overcomplete) set of exact GCS of Hj, sc \ and 
|0 > sc is the electron vacuum. 

In line with the general approach discussed in the previous section we can think of the 
i] k s as time dependent parameters, and construct the approximate trial time dependent 
wave-function of the full Hamiltonian Hn u b, \i>(t) >sc, as 

\m>sc=e^V\ri(t)> . (10) 
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We are now ready for evaluating the expectation value of Hsub over \if)(t) > sc , namely 
the semi-classical Hamiltonian 7"L r , which reads 



^-r 1 1 + rikTik 



+ VkVk 



VkVk 



+ VkVk 



k,l v 



VkVi 



VkVk)^ + vim) 



where i] k , fj k obey the Poisson-braket relations obtained from @, ih{v k ,V k } = (1 + VkVk) 2 - 
Instead of proceeding directly to the derivation of the canonical equations of motion, we 
notice that Hamiltonian ([11]) can be fruitfully rewritten in terms of the following semiclassical 
pseudospin variables 



S 



(k) 



1 VkVk - 1 

2 f] k i] k + 1 



z < 4>(t)\4 k} m) > 



s 



(fc) 



Vk 



; < 0(t)|jf |0(t)> 



1 + VkVk 



(12) 



and S_ = (S+')*, whose Poisson brackets recover for each k a su(2) k algebra. Explicitly 



2S 



(k) 



Moreover, one can define the related "mesoscopic" variables Sa^ = ^2 k£j \ <-'<> • 



S, 



(k) 



(13) 



with 



a = 3,+,—, A a = {k G A; e k = e a } denoting the mesoscopic (kinetic energy) levels. One 
can easily verify that the Sa a form a su(2) algebra like ([13]) (with k — > a), which we 
identify by su(2) a . Hamiltonian ([11]), when written in terms of Sa , reduces to a genuine 
one-dimensional problem, in that the index a (contrary to k) is strictly one- dimensional, 
numbering the different mesoscopic levels. Indeed 



H sc = 2j2(ta-fi)(Si a) + ^)+u 



S 3 + 



N 



+ |S 4 



(14) 



Here S a = ^2 a Sa^ = ^2 k Sa\ and u = U/N. The Casimir operators of both su(2) k and 

su(2) a algebras, I k = \S^\ 2 + \S$ | 2 (and the same definition for I a with k replaced by 

1 TV 2 

a) are conserved quantities for 7i sc . In view of the definitions (0), I k — and I a < — -, 
depending on the initial conditions. 

Noticeably, Hamiltonian (|T3|), like H and unlike H.[ sc \ commutes also with 5*3, i.e. the 
semi-classical variable corresponding to the total electron number operator. S3 is thus, as it 
should be, a conserved quantity, for which the relation holds 
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S 8 = y(n-1) . (15) 
In this sense we can therefore claim that the 1 — d Hamiltonian 7i sc obtained by means 

(sc) 

of the present semi-classical approach is a more accurate approximation of H than . 
In particular, in (|14|) the /c-modes are coupled dynamically through IS+I 2 , while in Ti!f c ^ 
they are not. This feature in turn keeps track in the present scheme of the non linearity 
of Huub, thus making Ti sc a good candidate for giving an approximate description of the 
physics of the Hubbard model in the whole phase-space. Of course, the results obtained by 

(sc) 

using instead H{ will be reproduced by the present approximation, as we shall see in the 
next section. 

From (p!3|) -([T4|) we can now derive the equations of motion for the mesoscopic variables 
Sa \ which read 



TthS { ± } = 6 a S { ± ] -2uS { 3 a) S± 



= u 



c c( a ) c c( a ) 



(16) 



Here 5 a = 2(e a —fi)+unN, where the constant factor unN—2/j, in S a is vanishing at half-filling 
(fj, = Tf, n = 1), and in any case does not affect the dynamics described by S± , apart from 
an over-all phase factor e ±: k( Un - 2 ^ T . Notice that of the three equations flT5"D only two are 
independent, whereas the third one is obtained from the Casimir constraint. For instance, 
one could use as independent variables S^ a \ which fixes also the absolute value of S±\ 
and the phase A a of = \S± \e ±tXa . This alternative representation of the pseudo-spin 
variables will also be considered, when useful, in the text. 

Let us emphasize that the true dynamical variables are of course the microscopic canon- 
ical variables Sl^, which satisfy the same equations of motion (0) with a replaced by k. 
Here we preferred to write them only for the mesoscopic variables because the Hamil- 
tonian 7i sc given by ([11]) was shown to be degenerate with respect to the inner dynamics 
of the mesoscopic variables. Moreover every solution we will be able to find for the 5±^'s 
holds straightforwardly also for the S± s, as 5k = S a . In fact, apart from this simple case, 
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every solution for the microscopic variables can be in principle worked out once we have 
found the mesoscopic solutions, and consequently S±, as in this case the equations flIS|) with 
a —>■ k reduce to a linear system with time-dependent coefficients. Interestingly, it is easily 
verified from ([TED that the scalar product of any two microscopic pseudospin vectors 
belonging to the same mesoscopic level is constant. This observation implies that in fact 
the time evolution of every microscopic vector in a given mesoscopic level is identical, the 
relative orientation of different S^ h \t) depending only on the initial conditions. 

According to the generalized TDVP approach introduced in previous section, and by 
means of (|3]), (|12D, (0), and (|16|), we are finally able to obtain the time-derivative of the 
time-dependent phase ip sc (t), 



A(fe) q (k) Mk) q (k) 

-H sc + ih22- 



1-2^ 



(k) 



2 N 2 
u < n — — 
4 



V l-2S< fc) 



-n sc + hj2^ + ^i k) )\ k , (17) 



where the last expression was explicitly written to make evident that a non-vanishing geo- 

(k) 

metric contribution to (p sc is expected whenever the phase of S± is not constant. 
Eq. (|17|) has some other relevant features which it is worth underlying: 

• it vanishes for vanishing u, as can be recognized from the second line form. In fact 
we know that if this is the case the wave-function given by ( |T0D becomes exact, and 
according to the discussion developed in previous section this implies that <^ sc (£) must 
reduce to the exact value given by which can be shown to be zero; 

,N 2 



it reduces to the constant un 2 —^- for S + = 0. Being S + related through ( |T2| ) to the 



semi-classical analog of the total pairing operator, it must be inferred that a non-linear 
time behavior of <p S c(t) is closely connected to the possible superconductivity of the 

state; 

contrary to both Ti sc and the equations of motion fll6|) , (p sc (t) can not be expressed 
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only in terms of the one- dimensional mesoscopic variables This means that it 

maintains memory of the inherent complexity of the original Hamiltonian, and gives 
information about the time-evolution of its wave-function which goes beyond that 



implicit in its semi-classical approximation (14), in particular depending on the inner 
dynamics of the mesoscopic levels a. 

IV. FIXED POINTS AND STATIONARY POINTS 

The first step in investigating the dynamical behavior of any nonlinear Hamiltonian 
system usually consists in finding its fixed points, that is those points in phase space where 
the equation of motions involve vanishing time derivatives of the dynamical variables. The 
stability analysis of such set of points leads to revealing their topological nature (by resorting, 
for example, to standard methods such as Routh-Hurwitz criterion) and, in conclusion, to 
structuring the phase space in regions where the dynamical behavior of the system exhibits 
well defined features^. 

A complete stability analysis is beyond the scope of present work. In fact, in this section 
we shall simply work out all the solutions to fixed point equations, in particular showing 
that indeed those among them which minimize the energy 7i sc give the same energy and 
the same self-consistency equation as the Hartree-Fock approximation, both for u < and 
for u > 0. Apart from that, the knowledge of the fixed points allows in principle to look for 
other solutions of (fl6|) by means of standard perturbative methods in their proximity. 

Minimum energy points are contained among the stationary points of 7i sc , which are 
easily shown to coincide with fixed points of ([16]) first by rewriting the equations of motions 
in terms of canonical variables (S3 , A a ), and by setting then A a = 0, S3 = 0. Since this is 
equivalent to = and = 0, eqs. ([L6D furnish the stationary point equations 

= 5 a S { ^ - 2uSi a) S + , (18) 

= si a) s_ -si a) s+ . 
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The case in which = for any a represents the simplest possible solution. We observe 

that S^ a \ thanks to the Casimir's constraint, can be chosen in a fully arbitrary way within the 

JV o (a) N a 

interval — < Sg < so that an enormous number of stationary points characterizes 
the mesoscopic pseudospin dynamics. 

It is important noticing that the solution S±^ = 0, when inserted in the equations for 
the microscopic variables, makes them immediately integrable, and the solution shows that 
in general a microscopic, inner dynamics for fc-pseudospins G A a can take place, according 
to 

S$\t)=S ( *\0)e a *t , (19) 

provided ^2 ke ^ a 5+^(0) = X]fceA a R e%Xk ^ — 0- Such a constraint in d = 2 is naturally obeyed 
by those configurations where the initial phase A&(0) is topologically nontrivial while R is 
independent of k. Indeed this is the case when Afc(0) -regarded as a function of k along 
the 1-d closed paths associated to each a-th mesoscopic level- undergoes a variation of 2np, 
with peN. For paths with energy e a ~ the number of modes iV a is great enough to allow 
e « A fe(o) t De twisted many times in a quasi continuous way. Let us underline that solution 
(Op, which -being consistent with S±^ = 0- corresponds to a stationary point of 7i SCl is not 
a fixed point of the microscopic dynamics when S±\o) ^ 0. 

The energy associated with the solution S+ = has the form E = — u(Ss + N/2) 2 + 
$0,(83 + F° r u > it is easy to check that an absolute minimum endowed with the 
energy 

Eit ) = -un 2 —-J2%\N a (20) 

a>F 

(a) N a 

is reached when = +(— ) — - for a > F (a < F) is imposed, F being that particular value 
of a for which 8p = 0, which implies fi = €f + U—, J2 a>F N a — J2 a<F N a = N(n — 1) —2S^ F \ 
and 5 a = 2(e a — ep). This absolute minimum corresponds to \S±^\ = for each k. Noticeably, 
the latter constraint does allow the ground state to still have a phase dynamics. In fact, 
on the one hand, by rewriting the equations of motions in terms of the canonical variables 
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A a , Sj^ introduced in previous section, it is straightforwardly verified that in the limit where 

(k) 

| S± I ^ uniformly the equation for the angle variables reduces to 

h\ k = 5 k -2uSi k) J2 C0S ( X k-*i) , (21) 

lek 

where, to second order in \S±^\, = ±| for k G A T . On the other hand for \S±\ = the 
phase is totally free, as one can check from equations fllq). Hence, by a continuity argument, 



we expect that also in this case the dynamics of Xk(t) evolves according to (pi). Equation 
(^T|) has already been investigated in a different context (see for instance^). In particular, it 
was shown0 that for XF-like models it allows for vortex-like excitations. Moreover, in the 
continuum limit it can be recognized as a Bernoulli-like equation, the latter being known 
to describe once more a vortex dynamics. Finally, let us observe that the solution of fl2T|) 
contains as a particular case the (topological) one discussed after eq. (|T9|), which requires 
S± 7^ 0, and reduces to it only in the (exact) non-interacting case, i.e. for u = 0. 

The absolute minumum ( |20| ) corresponds to the paramagnetic phase within the Hartree- 
Fock approximation, and gives the same ground state energy. Contrary to that approxima- 
tion, here it was possible to make evident a non-trivial dynamical behavior of the paramag- 
netic ground state. Such behavior implies in particular the appearance of a non-vanishing 
geometric phase in the ground state, as can be understood from the third of the equations 
(§). Let us recall that this should happen at any filling but half. We shall see in fact that 
at half-filling (n = 1) states built with antiferromagnetic order can provide lower energy 
for the corresponding semiclassical Hamiltonian, again in agreement with the Hartree-Fock 
approximation. 

Moreover, let us stress that stationary points characterized by S+ = -even when 
not identifying an absolute minimum- indeed can be shown to be local minima of the 
Hamiltonian when the geometric constraints represented by the Casimir's are taken into 
account. A simple first order expansion of 7i sc in the variables \S+ | 2 , where .Sg s are now 
expressed as = +(— )(/ — IS^I 2 ) -1 / 2 for a > F(a < F), shows that the variation 5H SC 
is positive provided u is positive and sufficiently small. In summary, we conclude that such 
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stationary points are minimum energy points for u > 0, possibly possess inner dynamics 
and topological structure, but do not involve superconductive situations, being S + = 0. 

The remaining set of fixed points - which are still solutions of equations (|T^) - can be 
fruitfully parametrized through the parameters I a and S + . Explicitly 



sA y^T4^p ' Si+a) ~ 2SaH \j 5l + Au"\S + \^ S+l ' (22) 

with s a — ±1. S + does not play the role of a free parameter, but it turns out to be 
constrained by the equation 



-*EW g+ £ w , ' (23) 



By substituting fl22|) in ([14]), and choosing the values of s a and J a which minimize 7i sc (s a = 1 
N 2 



and 7 a = ft is seen from ( p3|) that such solution exists only for u < 0, and corresponds 



to an energy 

J 2 
4 



yy2 i 

= M-7-n(2 - n) - - + 4m 2| 5+ |2 + | u || 5+ |2 (24) 



This result is once more in agreement with the ground state result for the superconducting 
regime within the Hartree-Fock approximation, as can be seen by identifying the varia- 
tional parameter A with the semiclassical pairing operator S + . In particular, the constraint 
equation (123]) coincides with the self-consistency equation for A. 



Let us notice that -as opposite to the mesoscopic fixed points S±* = case- here the 
insertion of the solutions fl22|) into the equations of motion for the variables S±^'s does not 
allow any microscopic dynamics, as the constraint S±^ ^ has to be satisfied. 



V. COLLECTIVE FREQUENCY DYNAMICS 

The dynamical system described by semiclassical eqs. (Tl6| ) is integrable in the special case 
when the pseudospin variables are supposed to be time-independent. The main effect 
of such an assumption is, in fact, of halving the number of the system degrees of freedom. 
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This can be easily seen by observing that |<S+^| cannot depend on time consistently with the 
fact that the Casimir's I a are constants of motion, so that only the phases of the pseudospin 
projection variables are allowed to depend on time. Further restrictions on the dynamics 
are due to the equations of motion which take the form 

-iKS^ = 5 a S { ? ) -2uS < f ) S + , (25) 

The first of eqs. (p5|) show how the system formally reduces to an ensemble of interacting 
oscillators with coupling constants 2uS^\ Moreover, together with the second, rewritten as 

w 

state that a unique, time-dependent phase — r characterizes the system dynamics. Namely, 

n 

for any a, 

S^(r) = V a e ia e^ T , (27) 

with V a , a, W G R. Both the constant phase a and V a are fixed by assigning the initial 
conditions (0). The linear character of eqs. (^) allows one to recast them in the matrix 
form 

{M{S} -Wt)-S+ = , (28) 

where the vector S + has components and the dynamical matrix M, whose elements 
can be obtained by system fl25|) , explicitly depends on the set S = {S^ }. The associated 
secular equation, which in turn provides the eigenvalue equation, 



det(M - Wl) 



q(b) 



h W-5 b 



U a (W-5 a )=0 , (29) 



is polynomial in W. The eigenvector components V a can now be expressed in terms of W, 
V = V a and as 

o,,c( a ) 



It should be noticed that in fact V is itself a function of the initial conditions and 



energy TC SC , through the relations ( 14|) and (0), which give V = ±[(TC SC + un 2 N 2 /4 — 

Moreover the eigenvalues fulfilling eq. (^9|) are obtained, after assigning the initial con- 
dition set {£3 }, by solving 



It turns out that the factor II a (W — S a ) in (|29| ) does not play a role unless S3 = for some 
a. When this is the case some of the eigenvalues coincide with the system proper frequency 
5 a . The total number of eigenvalues, corresponding to the number of different mesoscopic 
levels, is however kept constant. 

Both eigenvalues and eigenvectors can be easily obtained in an (approximate) explicit 
way when \u\/t is suitably small. Looking at the structure of (|3lD , it appears clear that the 
values of W close enough to 5 a are reasonably expected to fulfill it. In order to check this 



we first replace W with W a = e + 5 a in eq. (|31 ) which becomes 



5_ + 2u 22 



S b - S a - e 

then, by taking \e\ <C [ <5& — 8 a \, for any pair (a, b), one easily finds that e ~ —2uS^ thus 
obtaining 

W a ~5 a -2uSi a) . 

On the other hand, the condition |e| <C \5b — S a \ is satisfied if it holds in the less favourable 
case a = 0, b = ±1. Since |5 — 5i| = 2£ei ~ 8t7i 2 /N (for d = 2), then the condition on e 
becomes |e| <C 8tn 2 /N, which finally leads to 

|C/| < 8t7T 2 /iV a 

(a) 

for the greatest possible S% given by — . The present approximation scheme, based on 
considering W a ~ S a , is thus permitted for reasonably small values of \U\/t. In particular, if 
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\U\ <C 8t7r 2 /No all the eigenvalues can be obtained from this scheme, whereas if \U\ > 8t7r 2 
none of the W 's is well approximated by it. The a-th eigenvector associated with eigenvalue 
W a is readily obtained from eqs. (|30| ) and exhibits components V c (W a ) given by 



The above equations show how the a-th eigenvector S + (W a ) in this approximate case is 
characterized by the fact that only the component V& with b = a is strongly nonzero, being 
in fact V a ~ V. This implies that each eigenvector can be regarded as describing a supercon- 
ductive situation where the superconductive order parameter S + =< >= S+ is 
essentially given by S + ~ and the u-dependent contribution to the energy 7i sc is mainly 
given by the A;- modes with = e a . 

Some general observations are now in order. First we notice that the eigenvector problem 
is completely solved provided {S^ } and 7i sc -the quantities which, at this stage, describe the 
initial system configuration- have been assigned, and the eigenvalues W a have been worked 
out from (pip. No restriction constrains 7i sc and {S^} except for the filling condition fll5|) 
and the condition 7^0. The latter allows one to consider eigenvectors with arbitrarily 
small components V a but excludes the solutions characterized by S + = (V a = for 
any a) representing a subset of the solution of the fixed point equation (llq). The single- 
mode solution set is thus completely disjoint from such fixed point subset in the space of 
solutions of eqs. ( p5|) even if the former is dense around any element of the latter. On the 
contrary, the other fixed points of (p5|) , given by (^2|), are a (time- independent) subset of 
(p7[), corresponding to W = 0. 

Moreover we point out that the nonlinear nature of pseudospin dynamics survives our 
initial assumption = const because of the second of eqs. (0). In fact the linear 
system of coupled oscillators described by eqs. (|25| ) should have an arbitrary superposition 
of eigenvectors related to eq. (|28| ) as a general solution. This is no longer possible when eqs. 

are taken into account in that any superposition of single-mode solutions (eigenvectors) 
violates the request that pseudospins exhibit the same phase. 



V c (W a ) 



2uS£ 
5 C - 5, 



a 



V for c a 
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Furthermore we observe that S+ (r) can be easily obtained from eqs. (|25|) , where the 
term S+(t) is now playing the role of an external forcing term. Since S+\t) results to be 
proportional to up to a constant factor e l6> , then it follows that the equation for quantum 
phase ip S c(r) has form 

tUpsc(r) = W(S 3 + N/2)-nsc • (32) 

Hence, already in this simple integrable case within our approximation there is a non- 
vanishing contribution of the geometric phase (equal to WnN/2 times t), at any energy 
but the ground state. Such contribution should in principle be observable by appropriate 
experiments. 

Finally it is remarkable that these single-mode solutions, exhibiting some form of col- 
lective order through the unique time-dependent phase W-t, correspond to non-vanishing 
superconductive order parameter S+. 



VI. SLOW DYNAMICS VS FAST DYNAMICS 



A standard procedure for tackling many-body system dynamics consists in simplify- 
ing the equation of motions by separating fast degrees of freedom from slow degrees of 
freedomii'0. Such a procedure is profitable in that it leads the slow variable-system to 
become an autonomous system and sometimes reduces the complexity of its equations of 
motion. These features are, of course, appealing here because the slow variables dynamics is 
the one surviving at a macroscopic level (and thus it might be observable), while dynamics 
of high frequency degrees of freedom disappears on large time scales. 

For the pseudospin system a classification of pseudospins either as fast variables or as 
slow variables is naturally established by the fact that either S a ~ or 5 a ^ respectively. 
We recall that such reference parameter depends on n and is associated with the mesoscopic 
a-level of the ground state configuration where S^f 1 changes its sign. 
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The effect of such a distinction is made evident by performing the substitutions = 
exp(i5 a j:)ty a , which turn eqs. fllTf ) into the form 

-iMf a = -2uS { 3 a) y a -2uSi a) J2 ei{Sb ~ 5a) ^b , (33) 
ihS { 3 a) =u(v a J2 ne l(Sa - Sb ^) - c.c.) , 

b^a 

explicitly exhibiting dependence on the frequencies 5 a . Introducing the parameter 5 # as 
the frequency distinguishing slow frequencies (defined by \S a \ < <5*) from fast frequencies 
(defined by \5 a \ > <5*), it clearly results that those time-dependent oscillating terms of eqs. 
(|33l) where \8j,\ > <5* can be neglected on a time-scale greater than h/8*, since their rapid 
oscillations make their time-average vanishing. 

This fact has remarkable implications. In fact, upon denoting fast pseudospins variables 
and slow pseudospin variables by F± , and Q± , Q% respectively, we are now able of 
separating the dynamical equation set in two almost independent sub-sets, the first one of 
which describes short time-interval processes (t < h/8*), and reads 

-iHF® = 5 a Fi a) -2uFt\F + + Q + ) , (34) 



Fj a) (F_ + Q_)-c.c. 



while the second concerns long time processes (t > H/5*) involving the slow variables, and 
is given by 

- ihQf = 5 a Qf - 2uQf ] Q + , (35) 



thQ ( 3 a) = u{Q^Q 



c.c. , 



Here Q + = £ 'qS°, and F + = J2 a " F +^ where the 

prime and the double prime remind 
that a must range within selected intervals (\S a \ > 5* and \S a \ < 5* respectively). We notice 
that in eqs. fl34|) Q± can be regarded as time- independent terms (adiabatic approximation), 
since their evolution takes place on the time-scale of slow variables, whereas in eqs. (|35| ) fast 
variables are absent because of the effects of rapid oscillations discussed above. Also, when 
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t > h/5*, such oscillations makes F± negligibile with respect to Q± in terms like (Q± + F±) of 
eqs. (^5j), so that F- variable dynamics turns out to be driven by Q±. On the other hand the 
Q-system can be considered as an almost isolated system which exhibits the same features 
of the initial iV-pseudospin system except for the fact that now the pseudospin number is 
iV* = ^ a 'N a < N and the effective Hamiltonian is 

H Q = J2 2 (ta-riQ ( 3 ) +MN*/2 + Q 3 ) 2 + u\Q + \ 2 , (36) 

a 

where w* = U/N*. Two remarks are now in order. First we note that the long time dynamics 
is weakly influenced by those k-modes for which \5 a — 5f\ > so that the complexity of 
the dynamical behavior now issues from the Q-system, as manifestly suggested by the fact 
that the restricted Q-system has inherited the same structure the iV-pseudospins system. 
This is the main consequence of the adiabatic approach. Secondly, we recall that the density 
of states of the non-interacting system has in two dimensions a logarithmic divergence for 
which implies that the levels A a with e a ~ are the most populated ones. For 
situations where the value of n involves ep — (near half- filling) such a fact well matches 
with the first observation since it turns out that the Q-system, whose dynamics is complex, 
is also the sub-system involving the most part of fc-modes. 

From the above observations, one is led to restricting the number of interacting levels 
in order to work out the simplest yet still significant dynamics. The corresponding model 
turns to be a 3-level system, namely the pseudospin model where the Q-system is endowed 
with three levels. A simple calculation allows one to establish that the number of constants 
of motion is not sufficient to make the system integrable. In this sense 3-level dynamics still 
is far from being trivial, yet it is physically meaningful in several circumstances. 

At first, for example, one can take into account just the three inner most levels of the 
fc-space, i.e. those around e a = 0, which in the following we shall label by a = —1,0,-1-1. 
This is natural when investigating the low energy dynamics of 7i sc at half-filling with u > 0. 
In this case, in fact, it is reasonable to expect that increasing the energy from the ground 
state value of small amounts (recall that S$ = and = ±JV±i/2 with N\ = N-i) 
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makes interacting just the levels with the smallest energy, i.e. those with a = — 1, 0, +1. In 
view of the fact that |e a+1 — e a | <C 4i, expressing the almost continuous character of e a vs. 
a, one can replace both the upper level and the lower level of the 3-level model with two 
sheaves constituted by those levels with e a — e+i and e a ~ e_i respectively. This allows one 
to enlarge the number of modes partecipating to the dynamics as well as to treat situation 
where excited states are more than small perturbations. We recall however that this case is 
mainly pedagogic, as at half-filling the ground state of hamiltonian H has antiferromagnetic 
order. 

Such observations readily extend to those situations where 1. In these cases, in fact, 
the minimum energy configuration is not symmetric with respect to a = 0, but with respect 
to the level a = F, where <Sg change from negative to positive so as to minimize the energy. 
Then the 3-level construction must be referred to the new central level thus obtained. 

The further reduction to a 2-level scenario immediately makes the Q-system integrable. 
Again by replacing the two levels with two effective levels one can reasonably expect to still 
represent the main features of Q-dynamics, in particular when the energy is low enough to 
make interacting a limited number of levels situated around the level with a = F. 

VII. TWO-LEVEL DYNAMICS 

In the previous section we noted that the 3-level system is nonintegrable although it is 
the oversqueezed version of a multilevel system that was dramatically more complex. A 
thorough investigation of its dynamics, where the occurrence of a chaotic behavior indeed 
is expected due to its similitude with the dynamical model of refs. , requires a separate, 
extended analysis that will be pursued elsewhere. Nevertheless we shall start with the 3-level 
model equations, so as to make the approximations performed to achieve the 2-level scenario 
evident. 

Let us express the 3-level system equations in the form given by eqs. ( pop by renaming 
Sa (ct = +, — , 3) for v = F + 1, F, F — 1 by P a , Z a , and M a respectively, and by setting 
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5± = Sf±i, in order to simplify the notation and to recall the interpretation of the levels as 
level- sheaves. The equations then read 

- ihM + = <5_M+ - 2uM 3 Q + , ihM 3 = u(M + Q_ - M_Q+) ; (37) 
—ihZ + = —2uZ 3 Q + , ihZ 3 = u(Z + Q_ — Z_Q + ) ; (38) 
-ihP+ = 5+P+ - 2uP 3 Q + , ihP 3 = u(P+Q_ - P-Q+) . (39) 

For u > the 2-level model is readily obtained by freezing the ^-variables at the values 
Z 3 = = Z + . From the physical viewpoint, such an assumption is not particularly restrictive 
because it allows one to switch on dynamics through M-variables and P-variables starting 
from the ground state configuration. One should recall, in fact, that P + = M + = Z + = 
characterize the ground state, while the further condition Z 3 = can be implemented by 
suitably selecting n (see the discussion following (p0[)). The fact that the set S a is almost 
continuous still ensures the possibility of chosing n almost arbitrarily. 

For u < 0, instead, the fact that P + , M + , Z + ^ in connection with the minimum energy 
state, prevents the system from developing a dynamics in which Z + and Z 3 keep the their 
ground state values. Thereby the presence of the central level, even if as a nondynamical 
level reminescent of 3-level scenary, is prohibited and the central level must be embodied 
within one of the two sheaves, unless one is facing the nonintegrable version of pseudospin 
dynamics. At this point the 2-level scenario is restored and one can proceed to integrate the 
equations of motion. 

We construct now the solutions of 2-level dynamics by solving simultaneously the systems 
of eqs. (^) and (|39|). The main variable of the system is D 3 = P 3 — M 3 which will be 
shown to obey a nonlinear equation completely decoupled from the other variables. Indeed 
the knowledge of D 3 (t), together with the constant of motion Q 3 allows one to integrate 
such system^, which becomes linear with time-dependent coefficients. To work out the 
equation for D 3 one needs to exploit all the constants of motion. Explicitly, the energy 
H 2 = 5 + P 3 + 5_M 3 + u»(JV*/2 + Q 3 ) 2 + u\P + + M+| 2 , N* = N + + AL counting the active 
modes of 2-level dynamics, must be used to eliminate4 the variables P± and M± from the 

25 



equation 

h 2 D\ = 4u 2 |P+M_ - P_M + | 2 , (40) 

obtained via the second equation in (|37|) and in (|39|). This is done by exploiting first the 
identity |P+M_ - P_M+| 2 = 4|P+| 2 |M+| 2 - (|P + + M + | 2 - |P + | 2 - |M + | 2 ) 2 , which leads to 
rewrite the r.h.s. of (EOT) as 



.2| n 71 /r r> 71 /f 1 2 ic..2|d |2|7i/ |2 /ir, _ n „./lol2i|7i / r|2\ 



4xT P+M_ - P_M+r = 16u P+ M+ — 4 



ft _ aj D 3 _ u(\P+r + IM+ 



2 



(41) 



where h 2 = H 2 - 7 Q 3 - ^*(Qs + iV*/2) 2 , 7 = (8+ + <L)/2, and tr = (5+ - <L)/2. Then, by 
using the pseudospin Casimir's Im = M$ + \ M + \ 2 , I P = P 3 2 + |P+| 2 and the further identities 
P 3 2 + M| = (Ql + P 2 )/2, P 3 2 — M| = Q3P3, one reduces the P 3 -equation to the closed form 

^Dl = -2(h 2 -aD 3 ) 2 + 2u(h 2 -aD 3 )(2I-Ql-Dj)-2u 2 (I P -I M -Q 3 D 3 ) 2 , (42) 
with I = Ip + Im- 

Equation fl4"2]), which expresses the integrable character of the 2- level system, presents 
several interesting features. First of all it shows that the P 3 -dynamics is as complex as that 
of a 1-d potential problem. Indeed, upon introducing the potential 

U(D 3 ) = -2aDl + 26P 2 + 2cP> 3 + 2d 

where a = ua, b = a 2 + uh 2 + u 2 Ql, c = ua{21 — Q\) — 2ah 2 — 2u 2 Q 3 (Ip — Im) and d = 
(h 2 —ul) 2 + uh 2 Q 2 — Au 2 1 M I p , equation (f42|) simply becomes h 2 D 3 /2 = —U(D 3 ). Therefore 
the dynamical behavior of the 2-level system can be completely specified by identifying the 
regions where the cubic U(D 3 ) is negative and finding the value of the derivative of U(D 3 ) 
when P 3 approaches an inversion point. Such regions actually are identified by the compact 
interval in the potential well of U(D 3 ) whose extremes coincide with two of the three roots 
of the cubic equation U(D 3 ) = 0. The remaining semi-infinite interval where U(D 3 ) tends 
to —00 must be excluded in that D 3 there would assume infinitely large values, while its 
range is finite: — AL < D 3 < N + . 
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The two other points which play some role in characterizing the dynamics through IA(D$) 
are, of course, the potential stationary points 

6±(6 2 + 3ac) 1/2 j , 

furnished by dU/dD^ = 0. In fact, upon denoting the minimum and the maximum co- 
ordinates by D m and Dm respectively, it is possible to identify the fixed points as those 
configurations of the 2-level system with U(D m ) = so that the D 3 -interval reduces to the 
point D% = D m . It is easily checked that D m is consistent with the general result given 
by (f20|)-(f22|) for the exact minimum energy points. On the other hand, when the initial 
conditions imply U(D m ) < and U(Dm) > 0, then the system oscillates between two ex- 
treme states. From the physical point of view this property has the interesting consequence 
that the filling of each of the two mesoscopic levels near the Fermi surface (F ± 1) varies 
periodically with time, while the sum of the two fillings remains constant. 

On the other hand, it is worth noticing that such periodic behavior of the mesoscopic 
levels filling ceases to exist for appropriate initial conditions. In fact when U(Dm) = the 
system exhibits a laser-like effect, namely it tends, employing an infinitely long time, to 
an asymptotic stationary state in which (depending on the sign of u) a mesoscopic level is 
totally empty while the other is full. Moreover one should recall that each choice of the 
constants of motion h 2 , Q3, Im, and Ip embodied in a, b, c, and d selects a different cubic 
potential. Two solutions are associated with the same potential when they differ just for 
the choice of the initial position D^O). 

The nice feature of eq. (|42]) is that it can be reduced to the equation for the Weierstrass 
V function!!, which reads 



4V i -g 2 V-g 3 , 

Any solution of (H2j) can then be given in explicit form. A straigthforward calculation based 



on the substitution of D3 with D% = ±V + b/3a ( the plus (minus) corresponds to the case 
a > (a < 0)) turns ( [£2]) into the above equation for V where r = yfatjh and the standard 
coefficients g 2 , g% are identified as 
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<? 2 = ^(& 2 + 3ac) , g 3 = ±—(2b 3 + 9abc + 27a 2 d) . 



a 



Then, by exploiting the solution V in terms of Jacobi elliptic functions^ V(r) 
k 2 ) — 7 2 k 2 sn 2 (jt + a; k) where 





92 = ^7 4 (1 - k 2 + k 4 ) , g , = ±- 1 \l-2k 2 ){2-k 2 ){l + k 2 ) 



the explicit analitic expression of D 3 (t) is easily shown to 




(43) 



Such solution, as expected, shows that the dynamics is periodic around the point D m 
with period 

,K(k 2 ) 

since the elliptic sine fulfils the equation sn(x + 2K) = —sn(x), where K(k 2 ) is the elliptic 
integral of first kindil. As anticipated, a special case can be selected out when initial 
conditions allow the condition UiDu) — to occur. This condition states that one of the 
two W-roots confining the oscillations of D 3 inside the potential well, coincides with the 
maximum coordinate Dm- Reaching Dm thus requires an infinitely long time provided the 
motion starts exactly at the other inversion point, that is the remaining root of U{D 3 ) = 0. 
When this is the case k — > 1 so that sn(x, k) — ► th(x) and D 3 (t) describes indeed a transition 
for t — > oo. 

In order to esplicitly provide a situation where such transition happens we concisely 
examine the dynamics for u > when Q 3 = and Im = Ip- In this case the potential U 
manifestly exhibits its roots since it reduces to 



Hence the circumstance where U(D 3 ) = for D 3 = Dm is obtained by imposing the two 
roots of the quadratic factor in fl44]) to merge i.e. to tend to Dm- As a result one finds first 



U(D 3 ) = -2(h 2 - aD 3 ){h 2 - 2ul - oD 3 + uD\) 



(44) 
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the constraint 4w/i2 = a + 8u I on the energy, then that the limiting point of the transition 
is D 3 (oo) — ^- — D M when D 3 (0) = h 2 /a, provided the condition Ay/Tp < \a\/u < 2N + 
holds ensuring that — N + < D M < D 3 (0). 

Returning to (|4"3|), one can now easily evaluate the phase (f sc (t) given by expression (|I7|) in 



the simplest case in which S { + ] = S { + ] /N a for any k. Since this involves that S { 3 k) = ±S^ ] /N a , 
then pseudospins populating a level at most differ one from the other by the sign of Sf . 
We shall choose the same sign for the fc's of a given level so that the microscopic dynamics 
is just a copy of the mesoscopic dynamics. Here, we shall not consider the possibility of 
more structurated configurations of pseudospins, since they do not introduce any substantial 
novelty concerning the phase behavior. 

Finally, by using the above assumption, ip sc (t) can be written in the form 



-H sc + hY,(N a /2 + si a) )\ a . (45) 



The latter expression is particularly useful for the dynamics of weakly excited states since 
the variables S3 's are expected to undergo small variations in time, with respect to their 
ground state values. The levels whose pseudospins are fast should therefore contribute to 
ip sc principally through the phases A 's. On the other hand, after solving eqs. Q54] ) for the 
pseudospins labelled by a = F + 2 and a = F — 2 (namely pseudospins of the first two 
levels with fast variables) so as to have a 4-level system mimicking the real system, a simple 
calculation shows the time average of Xf±2 to be almost zero in the Q-time scale. Indeed 
the slow variables M3 and P3 and the phases of M + and P + provide the main contribution 
to if sc even when the fast variables are included. The expression of such contribution can 
be readily obtained by rewriting ([451 ) i n terms of M3, P3, and of the constants of motion. 
One thus finds 

The slow variation in time of its variables makes it a good candidate for experimental 
detection. 
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VIII. EQUATIONS OF MOTIONS IN THE ANTIFERROMAGNETIC PHASE 

When the possibility of a AF phase is considered, the natural order parameter which has 
to be non zero is0 

1 : X 

m = < Xy Gj ( n J<T ~ n id) >= AT < 5Z a k,T ak - G 't ~ 4a a k-G,l > , (47) 
jgA keA 

where G is a vector with all its components equal to ir. 

Also here, we look at the reduced Hartree-Fock Hamiltonian in the AF phase, in 
order to derive the dynamical algebra in which we shall subsequently construct our GCS's. 
It reads 



H} af) = ~t £ £4^ + |m]Te^K T - n u ) - U^(l - m 2 ) 

<ij> CT j 



2 v ">> ^' 4 

(«/) _ «•(«>/) 



and it can be rewritten in reciprocal space as H\ = X!keA_ a- His. a > wnere now 

H k,P = e k y^(wk l0 - - rik-G, CT ) + f^y (« kiT a k -G,T - « kii «k-G,| + h.c.) + j(m 2 - n 2 ) , 

and A_ is that half of A in which ek is negative, e.g. for d = 2 A_ = {k G A|ek < 0, or ek = 
OandO < h < n}. 

Hi can be recognized as an element of the dynamical algebra A a f = ■^■k,a '> with 

keA+ ,cr 

4? = [Kg = at G>a a kia ,K£ = [K$\\K™ = \{n^ - „ k , ff ) } ~ SM (2) k , CT , (48) 

and = K^J ± iK^J . In full analogy with the case treated in section III, we use 

as trial approximate time-dependent wave-function for studying the full Hamiltonian H in 
antiferromagnetic phase the time-dependent generalization of GCS's which can be built in 
A af , i.e. 

\m > af = e W*)|£( t ) >= eW*>n keA _ i(T (l +& i< £ wr )-iexp {^K^)\0 > af , (49) 

where |0 > a f= n ke A_ CT |cr > n kgi Y + CT |0 > (A + = A — A_), and the parameters £k,<x have to 
be thought of as time-dependent. 

30 



The semiclassical Hamiltonian TC a f, identified as the expectation value of H over 
\ijj(t) > a fi i can be fruitfully rewritten in terms of the semiclassical variables Ck,<r = a f< 
i/j(t)\K^\i/j(t) > a f, with C = X, Y, Z, which still satisfy a su(2) algebra (|T3|). Also here it 
is useful to introduce a one-dimensional index a instead of k, and to define the mesoscopic 
variables ( a)(7 = J2keA a ^ ne obtains 

H af = -2 eaZa^ + iuX^X^u— , (50) 

aeA_,cr 

With Ca = EaeA + Co,*- 

The semiclassical equations of motion for the Ca,o-'s, can be easily derived from (|50"D, (||), 
and read 



fcX acr 1t a Y a ]Cr , 



KY a ^ — — 2e a X aj(J — 4uX_crZ a ^ , 



KZ„ a = AuX_ a Y„ a 



(51) 



It is interesting to notice that the above equations do reduce to equations formally identical 
with those studied for the SC-paramagnetic phases for the special choice Y a = 0, X a = 
±X- a = 2^±> so that at least in this case the dynamics can be derived from that obtained 

±Y a ^ reduce to 



there. Moreover, let us notice that the choices X, 



±X. 



a half the number of equations (|5T|). One can easily verify that such choices minimize the 
value of H a f in the positive (— ) and negative (+) u regime respectively. 

From eqns. ([51]) we also obtain the time dependent phase charachteristic of the TDVP 
approach, 



Oaf = -Haf + Yl 



keA_ 



4«E 



ik,o-^"a-,(k) — -Xo-,(k)*k,<r 

1 — 2Zk,cr 

1 + 2Z Ka 

2Zk,cr 



keA_ 



(52) 



where 



N 2 



u- 



2 ek is a constant. As in the superconducting case, also in (|^) the 



keA_ 
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time-dependent part of ip a f is vanishing for vanishing u as well as for X a = 0, which is 
related to the vanishing of the antiferromagnetic order parameter X (X = Xj — X^). 

As in the case treated in the previous sections, also here we first look for the fixed points 
of eqs. (prrp. A first solution is of course the vanishing one, i.e. X a ^ a = F 0jCr = 0, and Z a>CT 



fixed by initial conditions. In particular, the configuration of Z a%a minimizing the energy 

N 2 

has energy E^j = —u h 2 e a N a , which is easily verified to coincide with that proper 



4 

a€A_ 



of the paramagnetic phase (^) in the positive u regime at half-filling. 

The remaining set of fixed points can be parametrized by the Casimir's X a a = X^ a +Z^ a , 
which again are conserved quantities. It is characterized by Y a>cr = 0, and 



X a , a -2s a , a u x l € 2 + ^2 X 2^ X -° > Za >°- eaSa '°\J el + 4u 2 X 2 _ a ' (53) 
where s aj(7 = ±1. The X CT 's have to satisfy the constraint equations X a = 



2u J2 a s ^\l t 2_^2 X 2 X_ a . The corresponding energy E af reads 



N / I 

Eaf = -U— + + + AUX ' Xi ■ 

In particular, the fixed points which minimize (El) are associated with the choices s a ^ — — 1, 
Za,a = and Xi = —sgn(u)Xi. In this case the constraint equations, apart from the 
solution X a = 0, reduce to one, i.e. 



\u 



V . , (55) 

V J el + 4u 2 X 2 



and the minimum energy is straightforwardly obtained from (|54|)-(|55|) as 



= -u^ - 2 £> fl) /e* + 4^X? + A\u\X 2 . (56) 

aeA- 

As expected, such energy corresponds to a non- vanishing antiferromagnetic order param- 
eter X = 2Xj only for sgn(u) = + (i.e. repulsive Coulomb interaction), whereas it gives 
X = for sgn(u) = — . In the first case, the energy coincides in fact with the one 
obtained within Hartee-Fock approximation, with X replaced by m satisfying the same self- 
consistency equation ([55]). On the contrary, in the attractive Coulomb interaction regime, 
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even though X = the energy Effl is lower than the one obtained within Hartree-Fock 
approximation, which would be precisely given by E^ , namely the energy corresponding to 
the trivial vanishing fixed-point. This is not surprising, in that while within the Hartree- 
Fock scheme the only parameter to be fixed self-consistently is m, here to all effect we have 
two related parameters, Xf and X[, which can separately be non-zero even when their dif- 
ference (i.e. X) is vanishing. Recalling that X a — a f< i>\ ) J ' n j,0-|'0 >o/> this latter case 
(X|- = Xi 7^ 0) can be recognized as a CDW phase. 

Notice that in the absolute minimum energy point for u < both the conditions which 
reduce the equations of motion ([y]) to those of the superconducting case fll6|) were fulfilled 
(see the discussion following eq. ([51]) ). A direct comparison with the result obtained for 
the negative u regime by means of the superconducting states (1241) shows that in fact at 



half-filling Esl ' = E^j . Hence we derived within TDVP at u < two degenerate wave- 
functions for the ground state, the superconducting and the charge- density- wave one. Indeed 
it is easily verified that the two wave- functions are orthogonal, and that the expectation value 
of the order operator of one phase, when taken over the wave-functions of the other phase, 
is identically vanishing. 

Now let us analyze the equation of motions fl5T|) away from the fixed points in some 
simple case. A first integrable case is obtained when the variable X a is kept constant. 
However this assumption is consistent only if ^ a e a F a](T = 0, and such condition in turn is 
satisfied only if Y a ^ is independent of time for e a ^ 0. Then the solution for each a ^ F 
reduces to (0), whereas for a — F it turns out to be given by 

Yp,<T — A a cos(a a t) + B a sm(a a t) , Zp >cr = A a cos(a a t) — B a sin(a CT t) , (57) 

and X F>a = A CT — ^ a _^ F A a (T , with a a = 4— X_ CT . Solution ([57]) survives in correspondence to 
stationary points of the Hamiltonian (when X a are chosen according to the self-consistency 
equations), as the system energy is not changed by the value of Zp^ a and Yf>. Such solution 
describes a periodic behavior of the mesoscopic Fermi level, holding even for the interacting 
ground state. In fact, due to the Casimir costraint the constants A a and B a turn out 
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to be related by the equation A\ + B\ = Xp — X Fa . The latter condition implies that 
A a = B a = for Xp >a = ±\/T~f, which is the case only for the absolute minimum point of 
the non-interacting case (see (|53"D). On the contrary, for any \X F>a \ < \fXp from fl5"T|) we 
obtain this oscillatory-like behavior of the solution at the Fermi surface. Such a behavior 
affects neither the order parameter, nor the energy, but it turns out to affect the phase (p a f(t) 
characteristic of the TDVP approach, by adding to the term linear in time a structured 
periodic time-dependent contribution given by 



E 

keA F 



tan 



B k ., T + v4 k)(T tan ( ^-t 



a„ 



1 - 4X2 



k*,<7 



(58) 



Here we used the same dynamics for the local and the mesoscopic pseudospin variables. In 
summary, we obtained also in the AF and CDW phases a non-trivial phase dynamics for 
the ground state. 

Apart from this simple case, more generally the system described by (|5l| ) has been inves- 
tigated in the case were the mesoscopic levels which have fast dynamics are one or twoHi. 
Already within such framework it appears to have very interesting chaotic properties. 



IX. CONCLUSIONS 

In the present paper we developed a consistent scheme for dealing with the dynamics of 
an itinerant interacting many-electron system described by the Hubbard Hamiltonian. Such 
scheme is based on TDVP procedure, and has been applied for describing the dynamics of 
the model by means of macroscopic wave-functions built in terms of GCS of the dynamical 
algebras which generate the Hartree-Fock solution in SC, AF, CDW, and paramagnetic 
phases. Already for these simple certain number of remarkable features related to 

the dynamical description rather than to the statistical-mechanical one was underlined. 

First of all, a geometric phase -a macroscopic quantity which in principle is observable- 
occurs for appropriate values of the physical parameters in the ground state as well as for 
some low-energy excited states . Such feature can not be identified by solving the eigenvalue 
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equation for the Hamiltonian (or related techniques, like the Bethe Ansatz approach), as it 
is a consequence of the phase of the eigenfunction, which in the eigenvalue equation is free. 
Even more noticeably, away from half-filling in the repulsive regime, it was shown that such 
macroscopic behavior of the Berry phase is originated from a vortex-like dynamics of the 
phases of the microscopic variables. Both these features could be due to the approximations 
implied by our scheme, hence a first interesting point which is left open to future work is 
to study exactly the dynamics of the stationary points of H, by solving the Schrodinger 
equation near them. This could be done by using the Glauber GCS's, which map exactly 
the quantum Hamiltonian into its semiclassical form, and studying the fixed points of the 
resulting equations of motionB. 

Other interesting dynamical properties of the system were stressed in the low-energy 
regime for some integrable cases. At half-filling, the ground state has been shown to exhibit 
an oscillatory-like behavior at the Fermi surface. Away from half-filling, for u > an 
analogous oscillating behavior for the mesoscopic density variable takes place near the Fermi 
surface. Such a feature is responsible for a non-trivial time dependence of the collective 
Berry phase. Again, this point should be further analyzed in different approximations. An 
alternative viewpoint could be furnished even by employing the same TDVP scheme starting 
from GCS's more realistic than the Hartree-Fock ones. For instance, in the U — > oo limit, a 
reliable basis is given by the Gutzwiller states^. 

One more solution obtained exactly within the present scheme and exhibiting interesting 
features is the single-mode solution, characterized by the non-vanishing of the supercon- 
ducting parameter, and by a unique time-dependent phase reflecting collective order. The 
possible relevance of this solution within the framework of superconductivity is related to 
the fact that it survives at u > 0, and at any energy but the ground state, with a volume in 
the space of solutions increasing with energy. 

All the above solutions, which are exact within the present approximation scheme, are 
interesting also in that they could represent good starting points for studying more ex- 
haustively the dynamics described by ( pf) in their neighbourhood, by means of standard 
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perturbative methods of classical dynamics. As a general conclusive observation let us notice 
that their validity beyond the present TDVP scheme could be tested by solving exactly the 
Schrodinger equation for Hu u b on small clusters of sites. Work is in progress along these 
lines. 
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